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Abstract 

In a gravitational virialized bound system built up of two components, one of which 
is embedded in the other, the Clausius' virial energy of one subcomponent is not, 
in general, equal to its total potential energy, as occurs in a single system without 
external forces. This is the main reason for the presence, in the case of two non- 
coinciding concentric spheroidal subsystems, of a minimum (in absolute value) in 
the Clausius' virial of the inner component B, when it assumes a special configu- 
ration characterized by a value of its semi- major axis we have named tidal radius. 
The physical meaning, connected with its appearance, is to introduce a scale length 
on the gravity field of the inner subsystem, which is induced from the outer one. Its 
relevance in the galaxy dynamics has been stressed by demonstrating that some of 
the main features of Fundamental Plane may follow as consequence of its existence. 
More physical insight into the dynamics of a two component system may be got by 
looking at the location of this scale length inside the plots of the potential energies 
of each subsystem and of the whole system and by also taking into account the 
trend of the anti-symmetric residual-energy, that is the difference between the tidal 
and the interaction-energy of each component. Some thermodynamical arguments 
related to the inner component are also added to prove as special is the tidal radius 
configuration. Moreover the role of the divergency at the center of the two sub- 
systems in obtaining this scale length is considered. For the sake of simplicity the 
analysis has been performed in the case of a frozen external component even if this 
constraint does not appear to be too relevant in order to preserve the main results. 
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1 Introduction 



In order to understand the dynamics of a two-component system wliicli a 
galaxy is, during its virial evolution, we need to explain what the trends of 
the Clausius' Virial and of the total potential energy of each subsystem, to- 
gether with that of the whole system, are, as soon as the inner (baryonic) 
component contracts inside the potential well of the outer dark matter halo. 
In spite of the limits of this model which, of course, is able to reproduce only 
some essential features of a real galaxy, its relative analytical simplicity, due 
to the use of tensor virial theorem, provides us with a powerful tool for under- 
standing the role of the physical parameters involved in the dynamic evolution 
of a real system. For the sake of simplicity we will assume that the outer com- 
ponent is at fixed size and shape without considering this constraint to be 
too essential in order to det ermine themain features of the dynamic evolution 
we are dealing with ( see, ( Seccol . 2001 ) , hereafter LSI). Indeed the masses 
of the two components are not equal, the outer one being about ten times 
the inner one. As a consequence, tidal influences between the subsystems are 
not symmetric; that acting from inner to the outer is actually weaker than 
the reverse (Caimmi & Secco 1993). Even if a contraction effect is induced 
from inner density distribution on the inner regions of outer halo, during the 
dynamical evolution, as already underlined by Barnes & White (1984) , that 
effect seems, indeed, does not cause a dramatic modification of the outer mass 
distribution, if there is a supernovae- driven outflow, according to some recent 
N-body simulations (Lia et al. 2000). On the other side, models in which this 
constraint has been changed with some less stringent additional conditions 
(Caimmi 1994), seem to prove the same conclusion. Then, in order to under- 
line the main effect of inducing a tidal scale length from the dark halo over the 
luminous component, we will neglect this possible modification. It seems to us 
that this may correspond to change from a given mass distribution to an other 
with a new exponent d (see, sect. 3) not too different from the previous one; 
that might be included in the different cases considered without cancelling the 
effect we are dealing with. 



The two subsystems will be modelled by two penetrating similar-strata spheroids 
with two power-law density profiles. The advantage of this model is to be able 
to handle all the energetics of a two-component system during its dynamic evo- 
lution in an analytical way. Nevertheless to obtain some physical meaningful 
relationships, we prefer to sacrifice the more correct, but very boring, exact an- 
alytical expressions and to use some suitable mathematical approximations. 
Through them it will be possible to describe the whole virial evolution by 
means of relatively simple equations which involve only the two exponents (6 
and d of sect. 3) which characterize the two subsystem power-law density dis- 
tributions (-B, the inner and D the outer, respectively), as soon as the common 
form factor F of sect. 3 is fixed (=2, in the spherical case). The appearance 
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of a non-monotonic trend for the inner Clausius' virial energy, the contrary 
of that which refers to the inner total potential energy, is the most important 
result for the dynamics of a two-component system during its virial phase. Its 
relevance has already been stressed by demonstrating how some of the main 
features of elliptical galaxy Fundamental Plane may follow as a consequence 
of it (LSI). Its occurence seems to be a general feature when one component 
is completely embedded in another. This has also been tested by consider- 
ing different models from those taken into account here a s such the case of 
heter ogeneous density profiles with different eccentricities ( Caimmi Sz Seccol . 
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Moreover, how much the extensions of the cores, at the center of the two sub- 
structures, have to be in order to obtain a maximum of the virial energy in 
the inner component, will be investigated when the luminous matter distri- 
bution becomes particularly steep. The present scheme may be extended to 
other two-component systems built up of baryonic + dark matter as, e.g., the 
galaxy clusters are. 



2 Virial and total potential energy 



By referring to Caimmi & Secco (1992) and the references therein, the Clau- 
sius' virial, or virial potential energy, Vu {u = B,D), appears with the kinetic 
energy, T^, in the scalar virial equations related to the two components, the 
inner, B, embedded in the outer, D, of a virialized system: 

2T, + K = 0. (1) 



It refers to the energy contribution due to the two active forces on the re- 
spective subsystem: the self-gravity and the tidal-gravity due to the other 
component. 

Generally speaking, it turns out to be: 

Vu = Qu + Vu,, (2) 



where 



r 3 

= J Pu / ^ Xr dXui (3) 
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$u and $t, being the gravitational potentials due to w-matter an v-matter 
distributions, respectively. 

In order to obtain the total potential energy of the w-component, we have to 
add at the self-potential energy the interaction energy Wuv, i.e.: 

{Epot)u = ^u + Wuv, (5) 

The difference between the total and the virial w-energy is: 

Vu — {Epotju = Quv, (6) 



Quv being the residual energy, which is anti-symmetric with respect to the 
exchange of one component with the other. The following equations also hold: 



Vuv ^^uv ~l~ Quvi (''') 
Wuv = Wuu, (8) 

the latter shows the symmetry of the interaction energy tensor. 

As consequence, the total potential energy tensor trace of the whole system, 
in the principal inertia axes frame of reference, is expressed by: 

Ept = {Epot)u + {Epot)v = nu + n^ + Wuv + Wuu = Vu + Vu. (9) 



It should be underlined that the residual energy may become zero, and then the 
total potential energy of each subsystem is equal to the Clausius' virial energy, 
only for some specific configurations which depend on the mass distributions 
of the two subsystems (see, sect. 7). 



3 Similar and similar-strata spheroids 



According to Caimmi ( 1993[ l (hereafter RC) and to ISeccol ()2000fl (hereafter 
LS), we consider, now, the case of two similar spheroids (that means, with the 



same axis ratio = = e (then the same form factor F = (2^^ -|- ^^^)) 
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with similar and coaxial strata mass distributions. We set the radial density 
profiles according to two power-laws as follows: 



PB~^, (10) 
1 , , 

PD^^; (11) 

and to avoid the central divergence we add an inner homogeneous core to both 
mass distributions. For the sake of simplicity we assume here that these two 
cores have the same extension (^c) in the a-dimensional coordinates (see, RC). 
Then, the traces of B and D self-potential energy tensors are, respectively: 



^B^-iyiiB ^F, (12) 

Ob 

QD^-i^nD^^F. (13) 

where is the mass and the semi-major axis of the considered subsystem. 
The coefficient vciu depends on the mass distribution via the functions h and 
gsu: from Table 1 of the Appendix A in RC. 

The traces of B and D tidal-potential energy tensors are: 



Vsu^-i^v^F, (14) 

as 

VuB^-l^VD^^F (15) 

an 

The symmetric interaction energy tensor trace is 

Wbd = -vw^^^F (16) 



The following expressions define the i/-coefficients of the previous energies: 



= -\ [{yB)M{yD)M\'^ mw''^\x), (17) 

^VD = -I [{yB)M{.yD)M]-^ — w*"*(a;), (18) 
8 mx 

^w^-^ [{yB)M{.yD)Mr^m[w'^\x)+w'-\x)l (19) 
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where x = as/ao e m = Md/Mb- The {i'u)m coefficients are defined in Tab. 
1 of RC. Due to the eqs.(17, 19) the following relation also holds: 

Wbd = i^[{yB)M{yD)M]-'m w'^\x)^^F + -Vbd (20) 
lb Ofi z 

The definitions of u''"*(x) and w^^*{x) and their approximations are given in 
the Appendix. 

It should be stressed that w^^* contains the expression of the Z)-mass fraction 
which has dynamical effect over the 5-component at every value of x (see 
Appendix) . 



4 Dynamical quantities 



In order to gain some physical meaningful relationships, we will adopt here 
some mathematical simplifications. It should be underlined that all the ap- 
proximations we consider in the following section and in the Appendix, refer 
to the case in which the two cores have the same a-dimensional extension 
and under the limitation that x > ^c- In ^ further paper we will relax these 
constraints. 

Under the approximations we have considered in the Appendix, eq.(14) be- 
comes: 

V,,.-^,G^F (21) 

Ob 



5-Clausius' virial may be defined as: 



Vb ~ GM^F 



1 , Md 1 
aB Mb an an 



<-y 



(22) 



or in normalized form: 

^^-GMIF- X ""^yMB^ ' ""-an ^ ^ 



By using the eq.(20) and the approximation (A. 14) (e ~ 0) we obtain: 

9r. X . X . ,GMbMd Cr^'^ 1 , , 
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and from eq.(7): 



Finally from the definition of {EpotjB (eq.(5)) and by the approximations of 
eqs.(24, 21), we have: 



{Epot)B^^B - -r[{'^B)M{'^D)M\ 



z(z — a) as CbD 
By normalization to the same factor of Vb (eq.(23)), it follows: 



or: 



where is: 



l.,_(3-c!)[5-(6 + <i)] , K 



M/B. = Wo. . (2- d)(3-6) + 2(2td) 



Moreover the residual energy of eq.(25) becomes: 

^ 1 / {?> - d)\h - ih + d)] (5-2d) , ^rf 



The Clausius' virial related to the D subsystem, eq.(2), is: 

Vd^^d + Vdb 

where fin and Vdb correspond to the eqs.(13, 15), respectively. In the limita- 
tions under which the approximation (A. 15) holds, the approximative tensor 
trace of Clausius' virial, related to D component, becomes: 

od od (2 — a)(3 — o) 2 — d^aD^ ' 
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By normalization at the usual factor ^MmZ obtain: 



(32) 



By using the definition of the total potential energy of D-component of eq. (5) 
and passing through the eq.(29), we obtain: 

/z? N 2 1 . {3 - d)[5 - {b + d)] , u{, 



The definition of total potential energy of the whole system, eq.(9), yields the 
following normalized and approximate expression: 

^"^"^ ^^"^ (2-d)(3-6) +(2^"""^ ^^^^ 



As a concluding remark of this section we wish open the question of which 
is the weight of the approximations on the trends of the energies if they are 
compared with those coming from the rigorous analytical expressions. This is 
remarkable (see the next section) when it refers to the non-monotonic charac- 
ter of the luminous Clausius' virial energy which turns out to be of extremle 
importance for the dynamics of the two-component system. Generally speak- 
ing, the approximations considered {b < 2) are as good as ^ (e.g., the 
maximum score for the main physical quantities of Tab. 1 is of 0.5% if = 0.01 
and b — 1.5), and that will constraint us to explore what happens on this limit. 



5 Luminous component energy trends 



We will analyze the trends of the two energies, the Clausius' virial of eq. (2) 

and the total potential energy of eq.(5) normahzed by the factor ^ , which 
correspond to the B component mass distribution (here assumed to be the 
luminous one). They are approximated by the eq.(23) and the eq.(28) respec- 
tively. The related rigorous functions are: 

t/^ ^ _^ _ ^ (35) 

X X 



and: 



{Epot)i 



X 



X 



(36) 
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Their corresponding trends are given in the Fig.l for the cases 2) and 4) of 
Table 1, by assuming = 0.01 and m = 8.5. From the comparison between the 
values of the main physical quantities of Tab.l, which are evaluated by using 
the rigorous trends shown in the figures, and those yielded by the approximate 
formulae, the conclusion is that the errors are less than 0.5%. This maximum 
error is reached only in the case 4). 

1 1 — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — I 








X 



Fig. 1. The energy trends of the i?-system as function of x = as /an, normalized 
at the factor {GMlF)/aD- The a-curve (solid line) represents the Clausius' virial 
energy, Vg, and the b-one (dotted line), the total potential energy, Eb-, (case 2) 
of Tab.l); a' and b' are the corresponding curves for the case 4) of Tab.l. For 
comparison the c-curve (long short-dashed line) is the self-energy trend of a single 
5-component in the case 2). The vertical marks correspond to the tidal radius 
configurations. 

It appears the relevant non-monotonic character of the Clausius' virial energy 
(curves a, a') in respect to the monotonic one of the total potential energy of 
the same component (curves b, b'). 

As already shown (see, LS and LSI), when we describe a gravitational struc- 
ture as in section 3, a configuration which minimizes the Clausius' virial of 
inner component exists, if the outer density distribution is described by a 
power law with an exponent less then 2. This special configuration of the 
embedded component is also analytically defined by the d'Alembert's Prin- 
ciple of virtual works. The condition of virial equilibrium means that, if the 
virial energy has a minimum (in absolute value), the total kinetic energy also 
exhibits a minimum. This implies, in turn, that if the kinetic energy of the 
inner bright component is dominated by the random velocity dispersion (as in 
ellipticals) the corresponding macroscopic pressure support has a minimum; 
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if, on the contrary, it is the rotational kinetic energy which dominates (as in 
the spirals), the rotation velocity distribution has to be that which is able 
to minimize the corresponding total ordered kinetic energy. The semi-major 
axis Qti which characterizes this configuration, is a scale length acquired by 
the gravitational field of the inner component, typically without any. It is 
expressed, by deriving the approximated eq.(23) in respect to x: 



1 



VOB 1 MbX'-' 



The result is the same as Hoerner ( 1958| ) found for the tidal radius of a spher- 



ical star cluster embedded in the galaxy tidal potential and represents a gen- 
eralization of it. This is the reason the dimension given by the eq.(37) is 
named as tidal radius. By assuming a^i = const., M^/Mb — 10 (on the basis 
of cosmological arguments) and by starting with the dynamical evolution at 
(^B = ctBo = cld {.x=1), we can see from eq.(23) that the ratio between the 
tidal-energy term {Vbd^ the second one) and the self-energy term (fi^, the 
first one) is: 

Ci = (^) - ^rn (38) 



For all the realistic values of the couple {b,d) we have considered (see, Tab.l 
and Fig.l in LS) the ratio (i is greater than one and then the tidal-energy 
term is dominant with respect to the self-energy one. As as decreases at 
ao = const., the tidal-energy decreases, in absolute value, because of the 
decrease in D-mass inside the i?-boundary, which is Mdx^~'^. Then the integral 
of force times the the position, Vbd, (21) decreases as meanwhile the self- 
energy increases in absolute value as x~^. The balance between the self- and 
tidal-energy is reached at about at and then, when x becomes less than Xt 
(= at/ an), the self-energy overwhelms the tidal-energy. The Ct ratio between 
the two energies at this minimum becomes: 

C, - (39) 



The approximated eq.(23) allows us to connect analytically the Xe value at 
which the two energies, self and tidal, are equal with Xj. Indeed, it turns to 
be: 

Xe~Xt(2-rf)3^ (40) 



li d = 1 the equipartition between the two energies properly occurs at Xt. 
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Neverthless the equation (37) is not the most general expression: the two 
components considered here are, indeed, always homothetic as as changes. In 
the homogeneous (6 = d = 0), non similar case, when the axis ratio depends 
weakly on a^, one approximative analytical expression of at exists (see, eq.(53) 
in LS). In order to look for a more general extension, as in the case of two 
different eccentricities with two heterogene ous density profiles, on ly a first 
order approach has been taken into account JCaimmi X. Seccnl .m; because 
in this case the tensor, (Kj„)ij {u,v = B,D; i,j = x,y,z), in an analytical 
form, is not yet available. 



It is interesting to consider where the tidal radius is located on the energies 
curves , in order to gain more insight into its physical meaning. Xt (the vertical 
mark on curves of Fig.l,) corresponds to the beginning of the self-energy 
domain in the Clausius' energy trend. Indeed, as soon as x < Xt both the virial 
and the total potential energy of the luminous component begin to increase (in 
absolute value) very steeply as the size decreases, with approximately the same 
trend of a single 5-component self-energy (curve c in Fig.l). At dimension 
greater than xt, the total potential energy curve of 5-subsystem exibits slopes 
smaller than the previous ones at a; < (see curve b of Fig.2), by changing 

its concavity a little bit before x = 1 (the zero of the c-curve(= ^^-^^5^) in 

Fig.2). If one looks at the derivative ^^'^^"'•'^ , curve b of Fig.2, it is manifest 
that at gives the scale length at which the slope of the B total potential energy 
changes. Indeed, if for x < Xt the slope is about the same of that of a single 
contracting component (curve b'), as soon as x > Xt the slope changes step 
by step in order to stabilize itself on one, which characterizes the tidal-energy 
domain, at about x ~ 3xt ~ 1. It should be noted that the flex point is included 
inside the interval Xt — 3xt (see Tab.l, even if in the case 4), strictly speaking, 
the mathematical flex point does not correspond to a physical configuration 
in our model because its coordinate is a little bit greater than 1). For x > Xt 
the rate at which the slope of the B total potential energy (curve c) varies is 
smaller than that at which the slope of same component without the potential 
well of dark matter (curve c') varies. 

In order to understand under which conditions the flex point does appear 
and where it is located, we will take into account the approximation given by 
eq.(28). Indeed, it allows us to handle the following analytical expressions. At 
flxed outer component, the first and second derivative of {Epot)B are, respec- 
tively: 



, ~ ^ + -uymx' (41) 
ax x^ 2 
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X 

Fig. 2. Detailed potential energy trend for the i?-system with and without the dark 
potential well, as function of x. The normalization factor is the same as in Fig.l. 
The a-curve (solid line) represents the total potential energy, {Epot)B in the case 
2) of Tab.l; the b (long dashed line) and c (dotted line) represent the first and the 
second derivative of the total potential energy, respectively, b' and c' are the same 
derivatives for the potential energy of the single B component. The vertical marks 
on the left correspond to the tidal radius configuration, the other, on the right, to 
the /Zex point of coordinate Xf (see, Tab.l). 

Then the flex point will occur at: 



It should be noted that the flex point may exist only if d < 1 and that it 
occurs at Xt times a factor / = (4^^)3-d (see, Tab.l) which depends only 
on the exponent d and which changes in a monotonic way from 2 to cxo as 
increases from to 1. 

We now turn back on the question, opened at the end of the previous section, 
related to the dependence of the non-monotonic character of the luminous 
Clausius' virial energy, and of the location of its maximum, on the approxi- 
mations used. As we have seen the approximated expressions of energies tend 
to the rigorous ones as soon as the 0. Then we have to explore what 

happens when the core dimensions in the two components decrease. In ac- 
cordance with our previous discussion on this point (see, sect. A. 2 in LS), we 
again conclude that, if the exponent b is not too high {b < 2) and the two 
power-law exponents satisfy the constraints given by (A. 8), the maximum on 
the Vb trend appears and its coordinate Xt changes very little as ^ 0. As 
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Table 1 

The most important physical parameters in the energy trends of figures. The com- 
mon curve parameters are: = 0.01, m = 8.5; xt = at/aD'-, Xf marks the ratio 
clb/o-d at which the variation of total potential energy of the S-component reaches 
its minimum (the flex, eq.43); x'^ marks the flex coordinate for the total potential 
energy of the whole system (eq.55). f = Xf /xt- 



cases b d xt f Xf Sxf x'j 



f 



1) 


0.0 


0.0 


0.389 


2.000 


0.778 


1.167 


0.590 


2) 


0.0 


0.5 


0.346 


2.702 


0.935 


1.039 


0.708 


3) 


0.5 


0.5 


0.361 


2.702 


0.975 


1.083 


0.739 


4) 


1.5 


0.5 


0.425 


2.703 


1.149 


1.274 


0.871 



soon as 2 6 < 2.5 the maximum also exists under the limitations (A. 8), but 
its location moves towards x — 1 with — ^ 0, by reaching a limit coordinate 
which is always less than 1. When 2.5 < 6 < 3 and again the hmitations 
(A. 8) are satisfied, the maximum appears if is not too small. Its coordinate 
Xt moves towards 1 as soon as decreases and then it disappears as ^ (in 
the case of, b=2.5, which happens at = 0.00038). These last two results are 
something new to add to the previous analysis. Then the following warning 
has to be pointed out: the non-monotonic character of Vb, may be lost ^ if the 
dimensions of the homogeneous cores become very small with 2.5 < b, even 
if the conditions (A. 8) are fulfilled. To have or not to have the divergency at 
the center of the two configurations, seems to become an essential feature in 
order to obtain or not an induced scale length on the inner component. This 
new message again requires, with the observations made at the beginning of 
the previous section, a further paper on this matter. 



^ It may be conserved from a mathemathical point of view, but the maximum 
corresponds to a non-physical configuration in our model, because it occurs at Xt > 
1. 
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6 Thermodynamical arguments 



To underline as special is the tidal radius configuration we will consider here 
some aspects of the thermodynamics related to the luminous component. As 
well kwown (Chandrasekhar 1939) when a single virialized component evolves 
through a quasi-static sequence of contractions, it has to lose an heat amount, 
AQ, which is one half of the variation of its self potential energy due to the 
contraction, in the meanwhile the other half goes to heat the structure. That 
means we need of a dissipation mechanism (e.g., gas cloud collisions) and of 
coohng processes. The consequence will be a local decreasing of the structure 
entropy and a corresponding increase of the universe entropy according to the 
II" thermodynamic principle (Secco 1999). Now, the question is: what hap- 
pens if we take into account the contraction process of the same component 
(B) if it is embedded inside an other non-dissipative one {D)7 In order to 
answer the question we have to re-consider the two fundamental equations: 1) 
the first thermodynamic Principle applied to B system and ii) the variation of 
the virial quantities of the same component during a quasi-static transition. 
In the 1° Principle equation it has to appear now the work done by both the 
forces on the system: the self gravity and the gravity which the dark mat- 
ter distribution exerts on it. The corresponding potential energy variation for 
a small contraction will be AVb- Combining the two equations and consid- 
ering that the internal energy is the total kinetic energy of the system, Tg, 
(Schwarzshild 1958), we obtain: 



where AQ is the amount of heat exchanged between the system and the uni- 
verse. Again the two requests that the system settles again in virial equilibrium 
and obeys to the 1° Principle yields the equipartition of the Clausius' virial 
energy variation that is one half of it has to be exchanged with the Univese, 
the other half has to contribute to change the "temperature" of the system 
which is proportional to the total kinetic energy. The result is formally the 
same found for the single component for which the Clausius' virial energy co- 
incides with the self potential energy, but the difference for the B component 
inside the D one is that the potential energy Vb now is a non-monotonic func- 
tion of X (Fig.l). The consequence is that the differential eq.(44) becomes the 
following: 



AQ = AVb/2 
ATb = -AVb/2 



(44) 
(45) 




(46) 



14 



where AQ is the exchanged amount of heat in normahzed units (sec. cq.23) 
and Ax the normahzed variation of the B structure semimajor axis. Moreover, 
by defining the entropy variation of the B system, in arbitrary units, as: 

= iAQ (47) 



The eqs.(47) and (23) yield: 



- p' m(2 - d)x^-'^ 
AS ^ ^ 



X 



+ Uymx"^'^ 



)Ax 



(48) 



where Tb = — |Vb, according to eq.(l). It is easy to see that the derivative of 
heat is stationary at x — Xt, and as consequence also the derivative of the B 
entropy. Then we have to wait a minimum or a maximum of the entropy of the 
luminous component at the tidal radius configuration. Indeed, the integration 
of eq.(48) yields always a non-monotonic trend for the function S{x) with a 
maximun at the tidal radius, as soon as the the maximun in Clausius' virial 
energy does exist. In the special case of Tab.l, b = 0;d — 0.5, we obtain: 

S(.) - = i„(12539) - fa( "339 + 1200. ^^^^ 



Its trend is shown in Fig. 3. 

The existence of this entropy maximum is able to stress the relevance of the 
Clausius' energy maximum. Indeed, according to Layzer (1976), the definition 
of thermodynamic information is : 

I = Smax ~ S (50) 



where Smax means the maximum value the entropy of the system may have as 
soon as the constraints on it, which fix the actual value of its entropy to S, are 
relaxed. Then, in order to increase its information a system has to decrease 
its entropy in respect to that of the universe. For a luminous component as 
B, which is embedded in an other D, a contraction decreases the entropy of 
B only if its dimension is smaller or equal to at- It means that only from 
this dimension downwards the component B may gain information during its 
evolution due to dissipation until to become more structured by changing its 
gas amount into stars (Secco 1999). Even if the argument has to make wider 
and deeper in a future work it seems to us it may become the ground for the 
interpretation of the oserved cut-off radius observed in many edge-on spiral 
galaxies (see, e.g., van der Kruit, 1979; Pohlen et al. 2000a, 2000b and Kregel 
et al., 2002). From a prehminary analysis the observed cut-off adius correlates 
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with the B mass as the tidal radius does (Guarise et al. 2001; Secco & Guarise 
2001). 



0.75: 
0.5- 
0.25: 
w - 
-0.25: 
-0.5- 
-0.75: 




X 



Fig. 3. Entropy trend of the i?-system, in arbitrary units, as function of x (case 2 
of Tab.l). The normalized factor is the same of the other figures. The maximum 
occurs at the tidal radius xt = 0.346. 



7 DM energy trends 



By looking at Fig. 4 which shows the energy trends of the D-component, the 
non-specularity of the virial energies behavior between the two components is 
manifest. Indeed, Vd is monotonic ( curve a), moreover {^EpotjD increases very 
slowly at decreasing x. 

It should be underlined that the two energies become equal as soon as the 
residual energies, given by the eq.(30), are zero. That occurs at the configura- 
tion characterized by the size ratio: 



(S-d) [5-{b + d)] 



(3 -b) (5 - 2d) 



(51) 



If the mass distributions of the two components are the same {b = d), the 
equality appears only when the two configurations are coinciding {xq = 1). In 
the cases of different mass distributions, it may be possible to find that the 
virial energy becomes equal to the potential energy only at one configuration 
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which satisfies the physical condition: < x < 1; in turn, that means : 

{2-d){b-d) < 

i.e., b < d (the case of the Fig.4), if we have taken into account the dynamics 
of two-component system at the presence of a tidal radius {d < 2). Then, the 
difference between the energies, at a; = 1, becomes: 



{Q 



BD)x=l 



iQ 



DB)x=l 



1 , (b-d) 



(52) 



All that is true not only for the D-component but, of course, also for the 
-B-one, due to the anti-symmetric nature of the residual terms. 
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Fig. 4. The trends of the D-system energies, as function of x (case 2) of Tab.l), 
normalized as in the previous figures. The a-curve (solid line) and the b-one (dotted 
line) represent the behavior of the Clausius' virial, Vd, and of the total potential 
energy {Epot)D-, respectively ; the c-curve (dash dotted -line) shows the trend of 
the residual energy Qdb = —Qbd- The two energies are equal at Xq = 0.9579 (the 
vertical mark) and the score at x = 1, is: {Qdb)x=i = 0.2361, eq.(52). 



8 The energies of the whole system 



In the last step we will investigate the potential energy trend of the whole 
system, which means the trend of the function (34) at decreasing x. It is shown 
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in Fig. 5. Again a flex point appears. Indeed, the flrst and second derivative 
turn out to be, respectively: 



dx^ 



dE. 



pt 



pt 



dx 



x" 



^ + u'ymx'-'' 
x^ 



(53) 
(54) 
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Fig. 5. The potential energy trend of the whole {B + Z?)-system, normalized to the 
usual factor. The a-curve (solid line) represents the total potential energy, Ept\ the 
b (dashed line) and c-curve (dotted line) represent the first and second derivative of 
the total potential energy, respectively. The curve parameters are those of the case 
2) in Tab.l. The flex point of eq.(55) is marked. 



Then the flex point will occur for: 



Xt 



(2 



2-d 1 

)3-c 



(55) 



It should be noted that again the flex point may exist only if < 1 and it is 
located a factor 1/23^ before the flex point of {Epot)B- 

The total potential energy trend of Fig. 5 for the whole system is also that of 
its total mechanical energy when it is in virial equilibrium. Indeed, according 
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to eq.(9), we have: 



Ept + Tb + T] 



'D 



Vb + Vd- \vb - -\vo = \Ept 



(56) 



The resuh is different for the single u-subsystem. From eqs.(5,6) we obtain: 



(-Smec)^ i^Epoi^u -\- Tu -|- Wu 



(57) 



The trends of the u-mechanical energies are given in Fig. 6. It is difficult to look 
for a physical meaning related to the small increase of the mechanical energy 
of the D component as x decreases (c curve). It is probably due to the artificial 
constraint, we have considered, that the dark matter system has a size and 
a mass distribution insensitive to the contraction of the baryonic component 
which occurs inside. As already pointed out by Barnes & White (1984), there 
is an iduced effect from B density distribution on the inner regions of dark 
halo which, even if non-dramatic, might change the coefficient uud, {t^d)m and 
then Uy. 
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Fig. 6. The mechanical energy trends of the u-subsystem, given by the eq.(57), 
normalized to the usual factor. The b (dotted line), c-curve (dash dotted- line) and 
a-curve (solid) represent the {Emec)B, {Emec)D and the residual energy Wbd, re- 
spectively. The curve parameters are those of the case 2) in Tab.l. 
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9 Discussion and Conclusion 



The energetics of a two-component virial system has been anahzed during its 
dynamic evolution due to the contraction of a inner (baryonic) component 
inside an outer of dark matter, which is assumed, for sake of simphcity, to 
be fixed. The results are also compared with those we obtain without the 
presence of the dark potential well. The analysis was done with the help of 
some simple, good approximate expressions for all the energies (in the cases 
of 6 < 2 and with = 0.01) which, once the form factor is assigned, are 
only functions of the b and d exponents in the two power-law density distri- 
butions. The different character between the energy trends of Clausius' virial 
and of the total potential energy for the baryonic (or bright component B), 
has been underlined and the physical main consequences has been derived for 
the minimum of the virial energy. The particular location of the special B size 
corresponding to the minimum {tidal radius) on the total potential energy 
curve and the discover of one flex point on it, has stressed the character of the 
tidal radius as a scale length induced by the dark matter halo component on 
the gravitational field of the inner baryonic one. This scale length marks the 
beginning of the two different regimes as follows: the self-energy regime, when 
X < Xt, and the tidal-energy regime at about Sxt- 

We claim as observable counter-parties of these theoretical reasons for the real 
existence of this induced scale length, at least the two foUowings: the first one 
is more direct: it is the cut-off radius observed in many edge-on spiral galaxies. 
Thermo dynamical arguments related to the presence of an entropy maximum 
at this scale length and some preliminary comparison between correlations 
from observables and theoretical expectations seem to prove the connection. 
The other one is more indirect but it has an important, observable conse- 
quence: if there is this length, we are able to reproduce some of the main 
features of the Fundamental Plane (FP) for the elliptical galaxies (LSI). From 
a general point of view, the two arguments are connected by the thermodynam- 
ics in the following way: the elliptical galaxies stay on the maximum entropy 
stage without a significant structure evolution, the cut-off spirals, practically, 
begins their evolution from this stage forwards. 

How the main physical results depend on the approximations used has also 
been considered. The main result is that: a too small homogeneous cores in 
both the two subsystems might cancel this induced length as soon as 6 > 2.5. 
Avoiding the divergency at the center of the two configurations, seems to 
become an essential condition in order to obtain the FP features. Here we are 
constrained by the restriction of handling together the two cores and then we 
need of a further paper on this matter in order to disentangle the probable 
different core-dimension play. 
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Nevertheless, from this point of view an appeahng apphcation to the galaxy 
clusters may be thought of. If we take into account the values of the exponent 
a, used by Girardi et al.(1995) in the power-laws which give the best fit for the 
galaxy distributions, at a distance r » {Rc= core radius), in the clusters, 
we can see that the corresponding range of our exponent b turns out to be 
approximately the critical one we have found in the section (5), 2.5 < 6 < 3, 
with ratios Rc/ Rvir (our ^c) which are big enough for obtaining the maximum 
virial energy at x < 1. 



The energies of the D-subsystem were also considered. Its total potential en- 
ergy is practically unaffected by the contraction of the -B-subsystem whereas 
the virial energy slowly increases, in absolute value, as x decreases. The con- 
ditions to obtain the equahty between these two energies was with the main 
result that only one configuration may enjoy this property. It corresponds to 
the initial one, when the two substructures fill the same volume, only if the 
two mass distributions are the same. 



Finally the trend of the potential energy of the whole system tells us that a 
configuration exists of minimum variation of the total potential energy (and 
then of the total mechanical energy), a flex point, inside its monotonic behav- 
ior. The conditions required in order to obtain either the flex on the B total 
potential energy , or on the total potential energy of the whole system were 
investigated and the result again is a constrain only on the dark matter dis- 
tribution (see the tilt of the FP). The mechanical energies either of the whole 
system or of both the subsystems were also taken into account. 



Even if we are aware that the next step we have to do is to investigate the 
energetics of a two-component system with or without central divergencies and 
characterized by two more realistic density distributions, e.g., of Hernquist 
(1990) kind for the inner component and of Zhao (1996) general kind for the 
dark matter halo (included the NFW profile (jNavarro et al. , 1993) ) , as already 
considered, e.g., by Caimmi & Marmo (2002), the present results seem to be 
very attractive for the future work. 
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A Appendix 



According to RC we define: 

w^"*(x) = L(ec,x) + M(ee,a;) (A.l) 
w'^\x) = G{^,,x)+H{^,,x) (A.2) 

dF 

L{^c,x)^ / Fn-r-^Cod^D (A.3) 
Fd-t^Cd dCo (A.4) 

^(ecx)^ / Fs-^eDdCD (A.5) 
{ d^D 

/dF 
FB-^iD din (A.6) 

The two functions Fb^b = ff'C-o) -Fd(Cd) in the ranges (0, ^c) and {^c, x), 
for different values of exponents b and are given in Tab. 3 of Appendix A in 
LS, with the assumption that (^c < f^)- 

In Table 1 of Appendix A (LSI), the corresponding functions Lij^^.x) and 
M(^c, x) are shown with the approximate analytical expressions of w"^^{x) we 
will use in the next sect.. In Table 4 of Appendix A (LS), the correspond- 
ing functions G{ic,x) and H{ic,x) are given with the approximate analytical 
expressions of w'^^^{x) which we will also need. We take into account the ap- 
proximation of w'^^^, towards which all the exact analytic expressions of it 
converge, by disregarding terms of higher order in respect to ^ different 
degree (see Table 4 in Appendix A of LS), which is: 

despite the b and d values, provided the following limitations hold: 

0<6< 3; 0<d< 2^{b + d) <5 (A.8) 



These limitations also ensure that the tidal and the self potential-energy tensor 
have the same sign, in all the cases we take into account. In order to high- 
light the D mass fraction, — Moi^)^"'^, which exerts dynamic effects 
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on B according to Newton's theorem (the equivalent of defined for the 
homogeneous case, see LS), we define the factor: 

q i:b+d 



In this way it follows that: 

(A.10) 

Mb an 



Moreover under the limitations given by the (A. 8), the following approxima- 
tion also holds (see LSI): 



w'^H^) ^ i^?-' ^ (^) + e (A 11) 

^ao' ~ {2-d)[5-{b + d)raD' (2 - d)(3 - 6) ^a^^ ^ ^ 

e^-A .^'^ ,.J —f5 (A.12) 



5 = 



(2 — d)2) ttB 
3 {aB\i> 



3 — 6 ^ttD' 



(^) ] (A.13) 



As we have already underlined (LSI) in case b = 0, 6 and then e of eq.(A.ll) 
become , whatever the value of d may be, and that the relevance of e, in 
respect to the other terms of the same equation, increases quickly, at fixed d, 
as soon as b increases. Moreover, at fixed b, the weight of e grows as d increases, 
even if less than before; at fixed b and d the contribution of e increases as the 
ratio ttB/aD becomes smaller. If we limit ourselves to choosing b < 1.5 and 
d ~ 0.5, we may disregard e in the approximation of ty*"* with a maximum 
error of 12% when b — 1.5 and aB/ao decreases to 0.4. That occurs if we 
choose = 0.1, but the error decreases as decreases. 

It should be noted that the approximation of w*"^* is different from that of 
yjext^ even if the core dimensions ^-re the same, because of the b value. 
Indeed, for 2 ~ 6 one needs to introduce in the approximated w^"'* the term S 
in which the ratio QB/ao has b as exponent. That does not appear in the w^^* 
approximation. 

By also taking into account the approximation of eq.(A.7) the following rela- 
tion holds, as soon as e is neghgible: 

^ A - w;-* (A 14) 

~ {2-d){3-byaD' ^™"(2-(i) ^ ' 
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By using eq.(A.14), (e ~ 0), we obtain the following approximation for the 
coefficient of eq.(18), as a function of z/y: 

.™ - (A.15) 

m 2 — 6 — b an 
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